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Abstract. The chemical evolution of galaxies that undergo an episode of massive and rapid accretion of metal-poor 
gas is investigated with models using both simplified and detailed nucleosynthesis recipes. The rapid decrease of 
the oxygen abundance during infall is followed by a slower evolution which leads back to the closed-box relation, 
thus forming a loop in the N/O-O/H diagram. For large excursions from the closed-box relation, the mass of the 
infalling material needs to be substantially larger than the gas remaining in the galaxy, and the accretion rate 
should be larger than the star formation rate. We apply this concept to the encounter of high velocity clouds 
with galaxies of various masses, finding that the observed properties of these clouds are indeed able to cause 
substantial effects not only in low mass galaxies, but also in the partial volumes in large massive galaxies that 
would be affected by the collision. Numerical models with detailed nucleosynthesis prescriptions are constructed. 
We assume star formation timescales and scaled yields that depend on the galactic mass, and which are adjusted to 
reproduce the average relations of gas fraction, oxygen abundance, and effective oxygen yield observed in irregular 
and spiral galaxies. The resulting excursions in the N/O-O/H diagram due to a single accretion event involving a 
high velocity cloud are found to be appreciable, which could thus provide a contribution to the large scatter in the 
N/O ratio found among irregular galaxies. Nonetheless, the N/O-O/H diagram remains an important indicator 
for stellar nucleosynthesis. 
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1. Introduction 

The nitrogen-to-oxygen abundance ratio in extragalactic 
Hn regions has long been found to increase with the oxy- 
gen abundance among spiral galaxies (Edmunds & Pagel 
1978, Vila-Costas & Edmunds 1993, Henry & Worthey 
1999 - Fig. Irregular galaxies are located in this plot 
near the lower end of the spiral sequence, and form a cloud 
with an appreciable amount of scatter around lg(N/0) = 
— 1.5 and over a range of oxygen abundances between 7 
and 8.5. 

Such a behaviour can be interpreted by galactic chem- 
ical evolution models as a signature of how the production 
of nitrogen in stars varies with metallicity. The closed-box 
Simple Model (neglecting the stellar lifetimes) predicts 
that the relation of an abundance ratio Zi/Z with the 
abundance of a primary element Z mirrors the metallicity 
dependence of the yield of element i, strictly independent 



of the star formation history and the specific recipe for 
star formation rate (SFR). This also implies that the rela- 
tion in the diagram can be interpreted as an evolutionary 
sequence from metal-poor to metal-rich systems. 

Thus, the N/O-O/H diagram has found a good use 
in infering the mode of nitrogen synthesis in stars. Since 
the ascending part of the diagram follows roughly a unity 
slope, it is taken to indicate secondary nitrogen produc- 
tion (i.e. with a yield proportional to the oxygen abun- 
dance). The horizontal part, where the irregular galaxies 
are found, then corresponds to the primary nitrogen com- 
ponent, i.e. with constant yield (e.g. van Zee 1998). Thus 
it appears that primary production dominates in the early 
phases of chemical evolution and secondary production in 
more evolved stages. 

This behaviour appears to be somewhat difficult to 
reconcile with what is known from stellar nucleosynthe- 
sis studies. While secondary nitrogen production occurs 
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in massive stars, an important contribution comes from 
intermediate-mass stars in the form of a predominantly 
primary production due to the Third Dredge-Up mixing 
freshly synthesised carbon into higher layers, where it can 
be turned by hot-bottom burning into nitrogen (e.g. van 
den Hoek & Groenewegen 1997). Thus, once the system 
is evolved enough for all intermediate-mass stars to con- 
tribute, the N/O ratio would remain rather constant in- 
stead of increasing with O/H. 

Several simplifications enter into such an interpreta- 
tion. One is that the lifetimes of nitrogen-producing inter- 
mediate mass stars are not negligibly short compared to 
the star formation. Such a delayed production makes the 
N/O ratio rise slower with time (and hence oxygen abun- 
dance) until it reaches the Simple Model relation, thus 
making the slope of the relation somewhat steeper during 
the early phases (e.g. Matteucci 1986). With this caveat, 
the N/O-O/H diagram is still an important instrument to 
probe stellar nucleosynthesis. 

The other, more important generalization of the 
Simple Model is the accretion of gas, whether by infall 
from outside the galaxy or outside the disk, or inflow 
within the disk of the galaxy. Koppen & Edmunds (1999) 
showed that as long as the ratio of accretion and SFRs 
decreases monotonically with time - which occurs in most 
chemical evolution models, e.g. with the infall rate falling 
off exponentially with time - the N/O-O/H relation re- 
mains within a factor of two close to the closed-box rela- 
tion. N/O ratios lower than that arc not possible in any 
kind of infall model. 

Detailed chemical evolution models - which take into 
account all the abovementioned details and processes - 
have confirmed that the N/O-O/H diagram still remains a 
very sensitive indicator of stellar nucleosynthesis. Henry et 
al. (2000) identified the horizontal part at lg(N/0) = -1.5 
as the result of primary nitrogen production in intermedi- 
ate mass stars. Because the lag times for nitrogen ejection 
are shorter than the long average star formation timescales 
in irregular galaxies, the N/O ratio reaches the plateau at 
lower oxygen abundances than in spiral galaxies. The rise 
of the N/O-ratio in more evolved systems is explained by 
the oxygen yield not being constant, but decreasing with 
metallicity. Such a behaviour would be expected from the 
evolution of massive stars whose substantial mass loss due 
to stellar winds increases with metallicity (Maeder 1992). 

There remains another potentially important aspect of 
accretion models: Koppen & Edmunds (1999) also showed 
that N/O ratios higher than obtained from the closed- 
box relation at the same O/H are possible, if the ratio of 
accretion and SFR increases so strongly and rapidly that 
the metallicity decreases in time. 

In recent years, it has become more apparent that in- 
tergalactic space is pervaded by neutral gas clouds. Among 
them are the isolated Compact High Velocity Clouds 
(CHVCs), defined to be kinematically not associated with 
a galaxy and not to have a stellar counterpart (Braun 
& Burton 1999, 2000, Kilborn ct al. 2000, de Heij et al. 
2002). Whether these clouds are left-over building blocks 



of galaxies (Blitz et al. 1999), remnants of galactic winds 
and fountains, tidal debris, or the product of ram-pressure 
stripping, their presence implies that they are likely to in- 
teract with galaxies. With typical sizes of 15 kpc and gas 
masses of several 10 7 M Q with clumps ten times smaller, 
they could affect the chemical evolution of the galaxy. 

One pertinent feature of the observed N/O-O/H dia- 
gram is the large amount of scatter, which is larger than 
the uncertainties of the observational material and the 
analysis methods (Henry & Worthey 1999). In particu- 
lar, one notes that while there seems to be a plateau at 
lg(N/0) = —1.5, there are more objects above this value 
than below. This tendency raises the question whether 
episodic infall the consequence of a collision with an in- 
tcrgalactic gas cloud - could contribute to the observed 
distribution of objects in the N/O-O/H plot. 

Over the past years much data on the gas content and 
the chemical composition of irregular galaxies of various 
masses has become available. This material will serve as 
a constraint to and a comparison with the chemical evo- 
lution models, in particular the relations between current 
gas mass and gas fraction, metallicity, and effective yield. 

In this paper, we first study the basic effects of rapid, 
massive infall on the secondary/primary relation (Section 
2) in order to show the basic behaviour and to identify the 
conditions under which strong deviations from the closed- 
box relation can arise. Section 3 presents the systematic 
relations among the observational data for a sample of gas- 
rich galaxies. Detailed numerical models for the chemical 
evolution of galaxies which undergo a recent infall episode 
are described in section 4 along with the results and the 
comparison with observational material. Discussion and 
Conclusions follow as sections 5 and 6. 

2. Effects of sudden gas infall on abundances 

The study by Koppen & Edmunds (1999) of the conse- 
quences of arbitrary infall of metal-poor gas on the chemi- 
cal signatures of galaxies showed that the determining fac- 
tor is the temporal behaviour of the ratio a(t) = A(t)/^{t) 
of the rates of mass accretion A and star formation 
If this quantity decreases monotonically with time, the 
galaxy evolves with a metal yield smaller than the true 
yield, but the abundance ratio of secondary and primary 
elements remains proportional to the primary abundance 
and is confined to a rather narrow region close to the 
prediction by the Simple Model Z p y s /y p > Z s /Z p > 

z P ys/(2y P )- 

On the other hand, abundance ratios larger than these 
limits Z s /Z p > Z p y s /y p can only be reached if the accre- 
tion is massive and rapid enough to cause a decrease of 
the primary metallicity. 

In the following we study more closely the basic be- 
haviour of the chemical evolution of galaxies that un- 
dergo such a massive infall of metal-poor gas. We con- 
sider models in which the finite delay between star for- 
mation and metal production is neglected. Oxygen is as- 
sumed to be produced primarily, with a yield yo = 0.0091, 
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Fig. 1. The abundance ratio N/0 as a function of oxygen 
abundance observed in spiral galaxies and irregular galax- 
ies after Henry & Worthey (1999), but with the symbol for 
Orion removed and that for the sun replaced by the large 
open circle. The dashed straight line is a Simple Model re- 
lation for purely secondary nitrogen production, while the 
top border of the shaded area is a relation for a mixture 
of primary and secondary production. 

to give results comparable to observed values. Nitrogen 
is taken to be produced purely in secondary mode with 
j/n = 0.095 • Zq. The SFR is supposed to be proportional 
to the gas mass, with a timescale of 5 Gyr. We note that 
this choice of the SFR does not determine the basic be- 
haviour in the N/O-O/H diagram, although it does affect 
the exact shape and timing of the evolutionary tracks. 

In Fig. |21 we show the evolution of nitrogen and oxygen 
of models which are closed-box, except for an infall of 
metal-poor gas starting at 5 Gyr age and lasting for 1 Gyr. 
It is evident that in order to have a significant deviation 
from the closed-box secondary-primary relation, the mass 
of the infalling gas must be substantially larger than the 
mass of the galaxy. 

Immediately before the infall starts, all models have 
a gas fraction of 0.55 and an oxygen abundance [O/H] = 
—0.3. The commencement of the infall raises the accretion 
ratio a from to about 10 and 1000 for M c i ou d/Af ga i = 1 
and 100, respectively. We note that in galaxies larger in 
size than a HVC, Af ga i will represent only that part of 
the galaxy which mixes with the infalling gas and forms 
stars. The oxygen abundance is strongly reduced due to 
the 'dilution' by the accreted metal-poor gas. As the gas 
in the galaxy accumulates, the star formation rate in- 
creases. Since the star formation timescale is longer than 
the duration of the accretion event, the enhancement of 
the SFR takes place appreciably after the infall. Thus we 
deal not with a proper starburst - in the sense of a re- 
duced star formation timescale - but rather with a galaxy 
which becomes more gas-rich and thus makes more stars. 
The oxygen synthesized by the massive stars limits the ex- 



Fig. 2. Tracks in the N/O-O/H-diagram of infall models 
with various mass ratios of the infalling gas and the galaxy. 
The infall event starts at the same age of 5 Gyr and lasts 
for 1 Gyr (marked by a dot). The shaded area indicates 
the region of the observations as in Fig. ^ 

cursion of the oxygen abundance, and it also brings down 
the N/O ratio. When the infall stops - marked by a dot 
on the tracks - the gas fraction is between 0.67 and 0.90 
for Mcioud/Mgai = 1 and 100, respectively. After this, the 
models continue to evolve like closed boxes, albeit with 
different initial conditions, and the tracks return to the 
Simple Model secondary-primary relationship. The over- 
all effect is that the model traces a 'loop' above the Simple 
Model relation in the N/O-O/H diagram. 

For the extreme model with an infall rate of 
100Mg a i/Gyr we follow in Fig. [3] the evolution along this 
loop. Because of the high infall rate the decrease in [O/H] 
takes place rapidly. Depending on the duration of the in- 
fall episode, the model returns within less than one star 
formation timescale to the Simple Model relation. If the 
infall episode lasts only a short time (e.g. 0.1 Gyr), the 
increased primary production of oxygen then leads to an 
evolution towards the lower right (similar to the zig-zag 
curves of Garnett's (1990) starburst model), followed by 
the (secondary) production of nitrogen leading towards 
the upper right, towards the Simple Model relation. If the 
infall lasts for a time comparable to the star formation 
timescale, the oxygen production may just compensate 
the continuing 'dilution' by the infalling gas. The oxygen 
abundance may remain nearly constant and thus one may 
obtain a rather large and deep loop leading again towards 
the close-box relation. 

The horizontal excursion of the loop depends essen- 
tially on the mass ratio of infall gas and the gas present in 
the galaxy. In Fig. 0] we study the influence of the time at 
which the accretion starts. As the gas fraction decreases 
from 0.80 (at 2 Gyr), 0.36 (at 8 Gyr) to 0.16 (at 14 Gyr), 
the loop makes a larger excursion when less gas remains 
available in the galaxy. 
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Fig. 3. The time evolution of models with an infall rate 
of 100 galaxy masses per Gyr, starting at an age of 5 
Gyr. The labels refer to the times after the start of the 
infall. The short-dashed curves show the evolution after 
accretion events of 0.1 and 0.3 Gyr. 



Fig. 5. Same as Fig.|21but with the proper nucleosynthesis 
recipe we will use in subsequent modeling. 



thesis prescription obey the characteristics predicted from 
the basic models: 




8.5 9 

ig(Q/H) 



Fig. 4. Evolutionary tracks of models with the same mass 
(10 galaxy masses) and duration of the accretion, but with 
different starting times (in Gyr, as labeled). All models are 
evolved up to 15 Gyr age. 



Finally, we present the evolution for the same model 
as in Fig. |3| but with the nucleosynthesis prescriptions we 
shall use in the further modeling (Section 0J: the produc- 
tion of nitrogen is taken from the stellar yields of van den 
Hoek & Groenewegen (1997). The SFR is assumed to be 
directly proportional to the gas mass with a timescale of 
5 Gyr. Because the primary production from intermedi- 
ate mass stars is important for nitrogen, the loops are less 
extended in the N/O ratio. This figure also demonstrates 
that the loops done with the full and detailed nucleosyn- 



— the N/O ratio does not reach below the limiting locus 
given by the closed-box model relation: this region in 
the diagram cannot be entered by evolution (Koppen 
& Edmunds, 1999). 

— during the accretion event, the N/O ratio does not 
become larger than the value just before the event. 
The abundance ratio cannot be changed by accretion of 
metal-poor gas, but only by the production of oxygen 
and subsequently of nitrogen. 

— the extent of the loop is governed by the mass ra- 
tio of the accreted gas and the gas left in the galaxy 
just before the event: If one neglects metal produc- 
tion during the event, infall of M c i ou( j of pure hy- 
drogen gas by a galaxy with a gas mass of Mi and 
metallicity Z\ leads to a reduced gas metallicity of 
Z 2 = Zi ■ M l /{M l + M cknid ). 

— the initial phase of a strong reduction of the oxygen 
abundance is rather short: it is determined by the infall 
rate which must be high if one wants a large excursion 
from the closed-box relation. Hence, it is more likely to 
observe accreting galaxies during the phase when they 
evolve back towards the closed-box relation. 

— short episodes of accretion lead to excursions mainly 
to lower oxygen abundances, while those comparable 
to the star formation timescale are necessary to make 
excursions towards lower N/O ratio: in short infall the 
gas is only 'diluted' whereas during longer episodes the 
resultant secondary (or delayed) production of nitro- 
gen follows only after the self-production of oxygen has 
started. 

We note that genuine starbursts - periods when the con- 
stant of the SFR is strongly enhanced - lead to a different 
characteristic behaviour, in particular that the N/O-ratio 
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reaches values below the closed-box relation (cf. the zig-zag 
tracks from Garnett, 1990) 

3. Observed relations among gas-rich galaxies 

From the literature we collect observational data on gas- 
rich galaxies. Data on abundances, Hi masses, gas frac- 
tions etc. are taken from Kobulnicky & Skillman (1996), 
Broeils & Rhee (1997), van Zee et al. (1997a, 1997b, 2001), 
Carignan & Purton (1998), Mateo (1998), Due et al. (1999, 
2001), Izotov & Thuan (1999), van Zee (2001), Kennicutt 
& Skillman (2001), Garnett (2002), Salzer et al. (2002), 
Young et al. (2003). Data for the same object from dif- 
ferent sources are combined, with preference given to the 
more recent data. 
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3.1. Mass-metallicity relation 

We show in Fig.Elthe relation between oxygen abundances 
at the effective radius and the Hi mass of irregular and 
spiral galaxies. There are 81 galaxies in this sample, and 
an overall regression can be obtained as 



12 + lg(0/H) = (4.93 ±0.37) 



(0.38 ±0.04) lgfe(l) 



which is almost identical to the relation obtained by 
Garnett (2002) 



12 + lg(0/H) = 4.31 + 0.43 lg(M H i/M ). 



(2) 



It is worth noting that spiral galaxies tend to have higher 
metallicities than this relation predicts for their masses. 
On the other hand, irregular galaxies with Hi masses 
larger than 1O 9 M have lower metallicities. 

3.2. Gas fractions and effective yields 

Garnett (2002) deduced gas fractions for 36 objects. For 
another 54 galaxies, we compute them from the ratio of 
Hi mass to blue luminosity and the extinction-corrected 
B-V colour index, using the recipe of van Zee (2001): 



j^gag _ Mm ^ g 1 10 2.84(B-V) -1.26 
Mstars L B 1.5 



(3) 



where the last two terms are based on the photometric 
properties of synthetic stellar populations. The gas frac- 
tions are clearly correlated with the Hi masses (Fig. I14|) . 
despite a large amount of scatter. It might be that this 
scatter is greater among dwarf galaxies. This behaviour 
is evident both from the galaxies with gas fractions from 
Garnett (2002) - marked by a cross - as well as in those 
whose gas fraction we computed from their Hi masses and 
blue luminosities. 

The Hi component in spiral and irregular galaxies is 
usually much larger in extent than the stellar disk. Thus 
the observed Hi mass may not be representative of the gas 
fraction in the galaxy. An example is Leo A where Young 
&: Lo (1996) find that emission from the inner galaxy is 



Fig. 6. Oxygen abundances and Hi masses of irregular 
galaxies of our sample. They are divided in three groups 
according to their Hi mass: small circles (Mm < 1O 8 M0), 
open circles (10 8 < M m < 10 9 M Q ), and filled cir- 
cles (Mhi > 10 9 Mq). Squares denote spiral galaxies. 
Wherever available, individual error bars are given. The 
dashed line indicates the relation by Garnett (2002). 
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Fig. 7. The gas fraction as a function of Hi mass in the 
galaxies in our sample. Symbols are the same as in Fig. 
El Objects with gas fractions from Garnett (2002) are 
marked with a cross 



about 20 percent of the whole emission. We show its po- 
sition in all figures for the entire Hi mass, but with an 
arrow pointing towards the loci if we took an Hi mass 
five times smaller. Fig. |Hl shows that this galaxy also has 
a very low ratio of dynamical to Hi mass, significantly 
smaller than most other galaxies. However, if one reduces 
the true galactic gas mass, this object becomes more sim- 
ilar to the others. 

With the knowledge of the gas fraction we can compute 
the effective yields, i.e. we compare the oxygen abundances 
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Fig. 8. The ratio of dynamical and Hi mass as a function 
of Hi mass. Symbols are as in Fig. El 
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Fig. 9. The effective oxygen yield relative to the solar oxy- 
gen abundance {Zo,q = 0.00959) as a function of the Hi 
mass of the galaxies in our sample. Symbols are as in Fig. 

El 

and gas fractions with the Simple Model relation. FigurcEl 
shows that despite an appreciable amount of scatter, there 
is a clear trend that more massive galaxies have higher 
effective yields. Spiral galaxies seem to be a natural ex- 
tension of the irregulars. This relation might indicate that 
the true yield is smaller in less massive galaxies, such as 
via a steeper IMF, but it could also reflect that the effec- 
tive yield is lower, because the stellar ejecta are less likely 
to be held back by the galaxies' gravitational potential. 

Figure shows the oxygen abundances plotted 
against the gas fractions, and compared to the Simple 
Model. It presents essentially the same information as 
Figs. an d El but for a smaller number of galaxies for 
which all three data are available. Among the irregular 
galaxies one notes a clear trend of more massive objects 



Fig. 10. The oxygen abundances and gas fractions (loga- 
rithmically on the lower axis, in percentages on the upper 
axis) compared to the Simple Model relations for solar 
yield (y = Zo,q = 0.00959, solid line) and one half, one 
quarter, and one tenth solar yields (dashed). Symbols are 
the same as in Fig. El 



having a higher effective yield, but the tendency for them 
to also have a lower gas fraction is much less evident than 
in Fig. Most probably this impression is due to the 
smaller number of objects and the large amount of scatter 
in all diagrams, some of which may well be genuine. The 
spiral galaxies are less gas-rich and have higher effective 
yields, and they might be interpreted as an extension of 
the irregular galaxies to higher masses. 

Finally, we show the N/O-O/H diagram for the galax- 
ies in our sample. For spiral galaxies we show the individ- 
ual Hn region data of van Zee et al. (1998b) in Fig. ^] 



3.3. Diameters 

Since the excursion from the closed-box relation in the 
N/O-O/H diagram is determined by the mass ratio of in- 
falling material to the gas still within the galaxy, we also 
need to consider what fraction of the galaxy will be af- 
fected by the infall of the gas cloud. According to Braun 
& Burton (2000), a clump in a high velocity cloud has a 
typical size of 1 kpc (a few arcmin at 1 Mpc distance). 
In a collison with a galaxy larger than this, only a frac- 
tion of the gas in the galaxy will participate and thus the 
chemical evolution is expected to be affected in this par- 
tial volume. To estimate this fraction, we need to know 
the size of the Hi component of the galaxy. 

Brocils & Rhec (1997) established a tight relation be- 
tween the diameter and the mass of the Hi disk in 108 
spiral and irregular galaxies: 

Dm = (Mhi/IO^Mq) 1 / 1 - 96 (4) 
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Fig. 11. Same as Fig.^ but for the galaxies in the present 
sample: spiral galaxies (small squares, individual Hn re- 
gions from van Zee et al. 1998b) and irregular galaxies 
(symbols as in Fig. H3 crosses: without Mm determina- 
tion). 



Fig. 12. The radial scale length of the stellar disk as a 
function of gas mass. The dashed line is the relation of 
Broeils & Rhee (1997) with the assumption of i?disk = 
D 25 /6. 



They also found a tight relation between the diameter of 
the optical disk and the Hi mass: 



D 



25 



(Mhi/IO^Mq) 1 / 1 ' 95 



(5) 



If we assume for the exponential scale length of stellar 
disk: 



Rdh 



D 25 /6 



(G) 



we get a very good match with the corresponding data 
of the objects in our sample, as shown in Fig. ^] Most 
of these data are from van Zee et al. (1997a, 2001) who 
find that the integrated luminosity profiles of the irregular 
systems are reasonably well fit by an exponential disk, 
and determine disk scale lengths by matching the surface 
brightness distributions in the R-band. 

van Zee et al. (1997b, 1998c) find that the ratio of the 
diameters seen in Hi and the optical -DH1/-D25 is about 2 
for 'normal' and low surface brightness dlrrs and between 
3 and 5 for BCDGs. There is no strong trend with Hi 
mass, as shown in Fig. ^] 

4. Models with detailed stellar nucleosynthesis 

We present calculations of the chemical evolution of galax- 
ies of a range of masses, using detailed stellar nucleosyn- 
thesis prescriptions and follow the evolution due to the ac- 
cretion of a clump of metal-free gas, taking into account 
that only the star-forming region is affected where the 
abundances are determined from the Hn region, so that 
in larger galaxies only a fraction of the galactic gas is in- 
volved. These models are constructed under the constraint 
of matching the various observational relations discussed 
above, and to show whether this process explains the scat- 
ter among the irregular galaxies in the N/O-O/H diagram. 
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Fig. 13. The ratio -DH1/-D25 of the diameter of the Hi 
distribution and the optical diameter. 



4.1. Assumptions of the underlying model 

For the prescriptions of stellar nucleosynthesis we follow 
Henry et al. (2000) and use the results of Maeder (1992) 
for massive stars and van den Hoek & Groenewegen (1997) 
for intermediate mass stars. We assume a Salpeter IMF 
with a stellar mass range from 0.1 to 100 M©. This gives 
a ratio lg(N/0) = —1.4 for the very low mass galaxies, 
slightly larger than the observed limit of —1.5. Although 
this offset is still well within the observational uncertain- 
ties and the accuracy of the nucleosynthesis results, it can 
be removed by slightly increasing the oxygen yield by a 
factor of 1.2. Thus the yield for oxygen (at solar metal- 
licity) is 0.026, about twice the solar oxygen abundance. 
Instead of introducing a galactic mass loss of some kind - 
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with its associated free parameters - we reduce both oxy- 
gen and nitrogen stellar yields by a common factor which 
depends on the galactic mass M ga io: 



V 

ynuci 



0.1 • M 8 - 4 if M 8 < 1 

0.1 • Af 8 01 if 1< M 8 < 31.6 

0.0355 -M? 4 otherwise 



with M 8 = A/ ga i, o /lO 8 M . The SFR 
1 ( 9{t) 



*(*) 



^SFR 



(g{t) + s(t)) 



g(t) + s(t) 



(7) 



(8) 



is assumed to depend on the instantaneous masses of gas 
g and stars s, thus on the total mass and the gas fraction. 
One advantage of this form is a well defined star formation 
timescale tsfr- In order to give a reasonable match of the 
observed dependence of gas fraction with gas mass, (Fig. 
01 we assume this relation 



tsfr = 20 Gyr ■ 



-o.i 



if Ms < 1 



M 8 otherwise 



(9) 



Here we do not assume the factor used by Henry et al. 
(2000) which enhances the SFR for higher metallicities. 
We also assume that the galaxy is formed by a slow con- 
tinuous infall with a rate of 



M{t) = M gali0 



exp(-t/rf) 
1 - exp(-T/r f ) 



(10) 



with a total evolution time of T = 15 Gyr and a timescale 
Tf = 4 Gyr, taken from the solar neighbourhood. This as- 
sumption is not crucial for our results: computation with 
a true closed-box model gives slightly larger oxygen abun- 
dances for the more massive galaxies. 

4.2. Modeling the accretion event 

From Braun & Burton (2000) we take that a typical 
CHVC has an extent of about 10 kpc (1 degree at 1 Mpc 
distance) with embedded clumps of about 1 kpc size (a 
few arcmin to 20'). The clumps have Hi column densi- 
ties of about 10 20 cm" 2 and thus masses of about 10 6 M Q , 
while the envelope (with 10 19 cm -2 ) is about 100 times 
more massive. For a typical relative velocity of the cloud 
of 100 km/s, this gives a duration of 10 to 100 Myr for an 
accretion of the clump and the entire cloud, respectively. 
It is much faster than the star formation timescale, esti- 
mated from the gas consumption, which is several Gyr. 
This implies that the deviation from the closed-box evo- 
lution cannot be a fully developed loop (cf. Figs. and [5J, 
but rather consists of an excursion to lower oxygen abun- 
dances, followed by the track evolving back towards the 
closed-box relation. 

As a simple approach the clump is assumed to affect 
only the region within r = 1 kpc of the galaxy's centre. For 
such a central collision, it is straightforward to estimate 
the mass fraction of the galaxy: 



M, 



gal, aft' 



M gal0 -(l-(l+r/R) 



-r/R 



(11) 



with the scale length R of the exponential Hi disk, which 
we take to be one sixth of the Hi diameter Dri- Combining 
this with the relation (Eq.EJ of Broeils & Rhee (1997) be- 
tween Hi diameter and mass, we get that while a gas clump 
would affect more than 80 % of the mass of galaxies below 
10 7 M©, it can influence only 14 % and 1.8 % in a galaxy 
with A/hi = 10 9 and 10 10 Af Q . This means that a gas mass 
of less than 10 8 Af Q is involved, and hence the accretion of 
a 10 8 Mq gas clump will make a significant impact on the 
chemical evolution of the galaxy's centre. For each galaxy 
we compute two models: one for the unaffected portion is 
done without an infall event, the other including the infall 
event. The stellar and gas masses of both parts are then 
added together. The metallicities are taken only from the 
affected part, since this will be what is measured in the 
current Hn regions. 

All galaxies evolve following the above presciptions, 
until the present epoch - 13 Gyr after the start of the 
evolution. At that time a metal-poor gas cloud is accreted 
over a duration of 100 Myrs with a constant rate. The 
mass of the cloud is randomly picked between 10 6 and 
lO 8 A/0 following a probability distribution proportional 
to M^ ud . We do not need to consider smaller masses be- 
cause of their negligeable effects. In all following figures 
we show the evolution of the models during the time in- 
terval between 13 and 15 Gyr after start of the evolution, 
i.e. what would be observable from now on. 



4.3. Results: gas fractions 

Comparison of the gas fractions (Figs. and I14JI show 
that the average relation among irregular and spiral galax- 
ies is reasonably well matched with our recipe (Eq. 0). 
The space between models with and without infall events 
would evidently be occupied by events involving less mas- 
sive clouds, which presumably occur more often. 

However, the models have far less scatter in the gas 
fraction than is evident from the observations. For masses 
less than about 10 8 ' 5 M Q the infall event causes some vari- 
ation of the gas fraction. But for higher masses the region 
affected by the accreting gas cloud is rather small, and 
thus the overall gas fraction is determined by the global 
SFR, which we had assumed to be a unique function of 
the galactic mass only. This assumption - taken for rea- 
sons of simplicity - does not appear to be generally valid. 
We take this as evidence of a genuine variation of the star 
formation timescale between individual galaxies, but we 
do not propose to model this here, because we wish to 
avoid introducing further free parameters. 

In Fig. one notes the presence of a few very low 
mass galaxies with low observed gas fraction: UGC 12613 
(type Im V, lg(M m /M ) = 6.61, lg(- ln(/ gas )) = 0.35), 
Pegasus (T = 10, 6.82, 0.033), and UGC 11755 (BCD/E, 
8.19, 0.46). They could be interpreted in terms of our mod- 
els as dwarf galaxies with a SFR higher than our recipe, 
having never undergone an accretion event, and thus con- 
tinuing to use up their original gaseous material. On the 
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Fig. 14. Gas fractions of the models as a function of their 
gas masses for times between 13 and 15 Gyr after the 
start of the evolution. The large plus-signs refer to models 
which do not undergo an accretion episode. The open cir- 
cles indicate models which are subject to an infall episode 
at 13 Gyr age. The symbols are evenly placed in time over 
the 2 Gyr time span. 



Fig. 15. Effective oxygen yields as a function of their in- 
stantaneous Hi - mass of the models during the 2 Gyr after 
the infall of a metal-poor gas cloud at 13 Gyr after start 
of the evolution. The horizontal dashed line indicate the 
solar yield. The scales of this figure correspond exactly to 
those of Fig. M 



other hand, there are the more massive irregulars with 
high gas fractions. The accretion event scenario might ac- 
count for them, if gas clouds of about 10 9 to 10 10 were 
involved. But they could also merely have a SFR lower 
than the other galaxies. 

4.4. Results: metallicity 

Our choice of the yield reduction factors (Eq. [7J| gives 
a satisfactory match of the effective yields of the mod- 
els (Fig. 1151) with what is deduced from the observations 
(Fig. EI) • The observed mass-to-metallicity relation is also 
reasonably well reproduced by the models (Fig. I17JI . 

The oxygen abundances and gas fractions, in relation 
to the Simple Models are shown in Fig. ^] Comparison 
with Fig. EH confirms that the choices for the dependen- 
cies of SFR and effective yields on galactic mass are quite 
reasonable. It also emphasizes that the irregular galaxies 
of moderate and high mass but large gas fractions (see 
Fig. 1 1 Of) are not well represented. While the two assumed 
relations could not model this, the episodic infall also can- 
not account for the presence of these galaxies. The low 
surface brightness "quiescent" irregular galaxy DDO 154, 
with its high gas fraction and high yield but otherwise 
with average properties (cf. Kennicutt & Skillman 2001) 
is a particularily striking example, as shown in Fig. 1101 

Our recipes for star formation and yield were chosen 
so that the resulting models give a reasonable match of 
the dependences of metallicity, gas fraction, and effective 
yield with the gas mass. The oxygen abundances for mas- 
sive galaxies thus are slightly higher than Garnett's (2002) 
mean relation, just as in the observations (Fig.HaJ). Among 
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Fig. 16. Oxygen abundances as a function of the gas frac- 
tion for our models. The full line indicates the solar yield, 
the dashed lines refer to twice, one half, one quarter, and 
one tenth solar yield. For comparison with observations, 
see Fig. UHl 



the dwarf galaxies, the models seem to indicate the pres- 
ence of many objects with very low metallicities. However, 
one must keep in mind that here the accretion of small 
clouds will dominate: events with clouds less massive than 
the minimum of 10 6 Af Q we show in the figures would fill 
up the space in the figures between the locus of the crosses 
and the model tracks with infall. As these events can be 
expected to be more frequent than collisions with larger 
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Fig. 17. The model oxygen abundances as a function of 
gas mass. The dashed line is the relation found by Garnett 
(2002). 

clouds, the emphasis in the diagram would be towards the 
less perturbed systems. 

4.5. Results: the N/O-O/H Diagram 

In the diagram of nitrogen-to-oxygen abundance ratio 
versus oxygen abundance (Fig. I18f) . excursions from the 
closed-box relation (indicated by the chain of plus-signs) 
are rather large but limited to a reduction in O/H in the 
horizontal part i.e. involving less massive galaxies. In the 
ascending part there are also quite strong effects due to 
more massive and hence metal-rich galaxies, which may 
populate the region above the horizontal part. 

If one looks at the evolutionary tracks in this dia- 
gram (Fig. I19|l . one notes that - by construction - low 
mass galaxies due to their longer star formation timescale 
evolve only until the horizontal part of the observational 
diagram, while more massive ones can reach up to the as- 
cending branch. For instance the model of 10 10 M© evolved 
(up to the present time) to a gas fraction of about 10% 
and an oxygen abundance somewhat larger than the sun, 
when the infall occurs. That a rather large loop is obtained 
has three reasons: firstly, a 10 8 M© cloud affects only 0.002 
of the galaxy's gas mass i.e. about 2 10 7 M Q which yields 
a reduction of the oxygen abundance by nearly 0.8 dex. 
Secondly, the star formation timescale is rather short - 
1.06 Gyr - and thus during the 100 Myr of the infall an 
appreciable quantity of oxygen is produced, and the N/O 
ratio can also be reduced (cf. Fig.|3J). Thirdly, the high star 
formation leaves a smaller gas mass in the galaxy at the 
start of infall. Note that during the infall episode the loop 
is passed through quickly, thus there is little chance to de- 
tect a galaxy in such a state. It is more likely to be found 
after the event when it evolves towards the closed-box re- 
lation. The enhanced density of symbols on these upward 
curved tracks makes them easily identified in Fig. El Since 



Fig. 18. The abundance ratio N/O as a function of oxygen 
abundance for the models. The symbols are as in Fig. 1141 
The crosses which denote models without infall are only 
barely visible and form the lower limit for the plotted cir- 
cles. The short-dashed curves envelope the maximum loop 
excursions for cloud masses of 10 7 and 10 8 M Q . The dashed 
straight line indicates the secondary nitrogen production 
track as in Fig. 2] 

the infall onto a less massive galaxy (e.g. 10 7 M©) affects 
the entire galaxy, the effects on abundances are also quite 
substantial, though not too evident in the diagram: the 
reduction of oxygen abundance is about 1 dex. But the 
evolutionary track after the event is at nearly constant 
N/O ratio, due to the slower SFR. This gives sufficient 
time for a stellar generation to eject all the delayed pri- 
mary nitrogen from intermediate mass stars along with 
the undelayed oxygen. 

As is evident from the tracks before the infall event, 
their slope becomes steeper in more massive galaxies. This 
is because the higher star formation rate makes the effect 
of the time-delay of the nitrogen production in intermedi- 
ate mass stars more pronounced. Thus, it is possible that 
the loop of a more massive galaxy may descent somewhat 
below the tracks of less massive galaxies. 

5. Summary and discussion 

Our model calculations demonstrate that it is possible 
to obtain a significantly broadened N/O-O/H relation by 
the galaxies evolving towards the closed-box relation after 
having undergone an episode of strong infall. 

One prediction of our scenario is that galaxies at high 
N/O ratio would tend to be more massive, as they origi- 
nate from galaxies which would have evolved into the as- 
cending branch, and the infall would have preserved the 
enhanced N/O ratio. Also, we expect that in the horizon- 
tal part and in the lower ascending part there is a mixture 
of less massive galaxies without an infall event which ex- 
perience a long evolution due to slow star formation and 
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lg(0/H) + 12 

Fig. 19. Complete evolutionary tracks in the N/O-O/H 
diagram: short dashed curves denote two models without 
infall events (galactic masses 10 8 5 and 10 10 ' 5 Mq); the oth- 
ers are models with a 10 8 Afo cloud falling on to a galaxy 
with masses of 1O 1O M (full line), and 10 9 , 10 s , and 10 7 
(with progressively shorter dashes). 

of more massive ones having undergone an infall episode. 
Comparison with Fig. 1111 shows a substantial amount of 
mixture, but the observational uncertainties in the abun- 
dances and the lack of Hi mass determinations do not 
permit to clearly confirm or disprove the expected trends. 

5.1. The N/O ratio 

The computed N/O ratios as a function of galactic Hi 
mass, presented in Fig. 1201 are nearly constant for small 
masses, since the loops due to infall events take place in 
the fiat part of the N/O-O/H diagram (cf. Fig. 0. Only 
for high mass models does the ratio increase by about 0.5 
dex. This is because the massive objects would cover the 
ascending branch. Here, infall episodes may lower the N/O 
ratio, but only by about 0.2 dex. 

Among the observations (Fig. I21fl we note that apart 
from II Zw 40 and the Pegasus dwarf, the irregular galax- 
ies show no clear trend with Hi mass, but seem to scatter 
around a constant value. In this figure none of the objects 
marked with crosses in Fig. II II is present. The Hn region 
data from spiral galaxies show the range of N/O ratios 
met in these objects. On average, the values are slightly 
higher than among the irregular galaxies. 

For II Zw 40, Kobulnicky & Skillman (1996) give 
lg(N/0) = -1.05 while Izotov & Thuan (1999) give -1.53 
which is more like the other objects. The value for the 
Pegasus dwarf of —1.24 is found by Skillman et al. (1997) 
to be one sigma above the values of the other irregular 
galaxies. For this exceptional object, the authors discuss 
several possible explanations, including a different IMF 
and time delayed nitrogen production. This would remain 
valid in the context of episodic infall. We shall not consider 




7 8 9 10 

ig(M gas [M SUN ]) 



Fig. 20. The N/O abundance ratio as a function of gas 
mass, as computed from the models. 




_2 L , i , , , , i , , , , i , , , , I , ._ 

7 8 9 10 

lg(M HI /M SUN ) 

Fig. 21. The observed N/O abundance ratio as a function 
of galactic Hi mass. The data from spiral galaxies (small 
squares) pertain to individual Hn regions, while the Hi 
mass of the parent galaxy is used. The other symbols are 
the same as in Fig. 

these two cases. We also note that our models with their 
single dependence of parameters on galactic mass cannot 
account for the appreciable scatter in the diagram. 

5.2. Models with additional starbursts 

To explore the additional effects due to the presence of 
genuine starbursts, we compute the same series of models 
with a bursty star formation history. We apply the fol- 
lowing ad-hoc recipe: for all galaxies less than 10 10 Af© 
the star formation rate is 10 times the continuous rate 
used before, for a duration of 200 Myr beginning every full 
Gyr, but during the 800 Myrs in between the rate is one 
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lg(0/H) + 12 

Fig. 22. The N/O-O/H diagram for models with addi- 
tional genuine star bursts. 
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Fig. 23. Evolutionary tracks in the N/O-O/H diagram for 
a model with 10 9 Af© mass and with genuine star bursts. 

hundred times less the continuous rate. This increases the 
scatter in the N/O-O/H diagram (Tig. 1221) somewhat. One 
notes that the closed-box relation is no longer as strictly 
obeyed as a limit; as expected, several points lie below the 
relation. Likewise, in all the other diagrams the amount 
of scatter among the dwarf galaxy models is further in- 
creased, although this does not significantly change the 
relations seen in Figs. El t° El The evolutionary tracks 
(Fig. I23|) before the infall event show the typical Garnett 
zig-zag curves, although somewhat modified because of 
the continuous infall. The accretion event gives a large 
loop, with a variation of about 0.1 dex in the N/O ratio. 

6. Conclusions 

We examine the chemical evolution following the rapid 
infall of metal-poor gas into a galaxy. First, we consider 



models with a simplified nucleosynthesis recipe and ex- 
plore the influence of the various parameters, and under 
which circumstances one achieves large deviations from 
the closed-box Simple Model. We find that 

— excursions in the N/O-O/H diagram are in the form 
of loops but only towards the 'upper-left' side of the 
closed-box relation, that is towards smaller oxygen 
abundances and larger N/O ratios. 

— during the infall, the oxygen abundance is decreased 
due to dilution of the galactic gas, followed after the 
infall by the evolution towards the closed-box relation. 

— in order to have large excursions, the mass of the in- 
fall gas must be much larger than the mass of the gas 
present in the galaxy (or in the star-forming region of 
large galaxies). 

— also, the infall rate must be greater than the SFR. 

— in order to form large loops, the duration of the event 
must be of the order of the star formation timescale. 

These results, which may also be obtained by analytical 
solutions, are also found from numerical models with a 
detailed nucleosynthesis recipe and taking into account 
the finite stellar lifetimes. 

We then examine the possible effects of HVCs collid- 
ing with gas-rich galaxies on the chemical evolution of the 
galaxy. From the properties attributed to HVCs we esti- 
mate that 

— a passage would occur for 10 to 100 Myr which is suit- 
ably short compared to the star formation, 

— the gas masses of 10 6 to lO 9 ^/© would be sufficient to 
cause a large influence on dwarf galaxies, 

— in massive (and large) galaxies the HVCs size of 1 to 
10 kpc would restrict any observable effects to a partial 
volume of the galaxy. 

This implies that such a collision could be expected to 
leave observable marks in the chemical properties of the 
galaxies, different from what a closed-box evolution pre- 
dicts. Based on the experiences of the simplified models - 
but confirmed with the detailed models, the above char- 
acteristics also imply that 

— full loops in the N/O-O/H diagram are not realized, 
because the infall duration is short, 

— large loops or excursions will be rare, since they require 
the encounter with a massive cloud component which 
most probably is a rare species, 

— because of the short infall duration, galaxies are most 
likely to be observed during the slow evolution back 
towards the closed-box relation. If infall triggers en- 
hanced star formation, the galaxies might be very 
bright at that early stage. 

To model the global properties of galaxies of various 
masses, we assumed relations of star formation timescale 
and yield on galactic mass which we adjusted to repro- 
duce the observed average relations of gas fraction, oxygen 
abundance, and effective yield with Hi mass. Evidently, 
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our models cannot account for the appreciable dispersion 
of the observed galactic properties, because e.g. infall hap- 
pens in reality during all stages of evolution. Following the 
chemical evolution for the 2 Gyrs after a single encounter 
with a HVC, we find that among galaxies less massive 
than about 10 8 M Q the deviations from models without 
infall events are substantial enough to provide some con- 
tribution to the observed scatter of galactic properties. 

In the N/O-O/H diagram, the galaxies are mainly 
caught on the slow way back to the closed-box relation. 
Bright starburst galaxies might be observable at earlier 
stages, such as near the largest decrease in oxygen abun- 
dance. All this adds a quite significant dispersion to- 
wards the 'upper-left' side. Thus the dispersion of irreg- 
ular galaxies towards N/O ratios higher than the plateau 
at lg(N/0) = —1.5 might well be objects that had evolved 
until higher metallicities than they now have, but which 
suffered an infall episode with a HVC. 

On the other hand, we note that the deviations remain 
somewhat limited; galaxies are likely to be observed in the 
vicinity of the closed-box relation. Thus, the N/O-O/H 
plot retains its sensitivity to the way the two elements 
are produced and incorporated in the interstellar gas. The 
diagram remains a very useful tool to infer details about 
stellar nucleosynthesis. 
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